Rung-singlet phase of the S = ~ two-leg spin-ladder with four-spin cyclic exchange 
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Using continuous unitary transformations (CUT) we calculate the one-triplet gap for the antiferro- 
magnetic S = \ two-leg spin ladder with additional four-spin exchange interactions in a high order 
series expansion about the limit of isolated rungs. By applying a novel extrapolation technique we 
calculate the transition line between the rung-singlet phase and a spontaneously dimerized phase 
with dimers on the legs. Using this efficient extrapolation technique we are able to analyze the 
crossover from strong rung coupling to weakly coupled chains. 
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After the discovery of high-T c superconductivity in 
1986, low dimensional quantum antiferromagnetism has 
attracted much attention in condensed matter physics. 
Recently it has become clear that the minimum magnetic 
model for cuprate systems has to include four-spin ex- 
change terms besides the usual nearest neighbor Heisen- 
berg exchange interaction 1 " 10 . One important subclass 
of such models are the two-leg ladder systems. The 
nearest-neighbor Heisenberg model on the two-leg lad- 
der without four-spin interaction is a gapped spin liq- 
uid. This system is in the rung-singlet phase and first 
excitations are triplets 11-14 . In the limit of zero rung 
coupling there are two isolated gapless spin chains. In- 
cluding four-spin exchange interactions several new quan- 
tum phases occur 15-17 . Possible phases include a spon- 
taneously dimerized phase where the dimers are located 
in a meander-like structure on the legs, scalar and vec- 
tor chirality phases, a region of dominant collincar spin 
and a ferromagnetic phase 16 . However, real two-leg lad- 
der cuprate systems are always in the rung-singlet phase 
but relatively close to the quantum phase transition to 
the spontaneously dimerized phase 4 ' 9 . Therefore, it is in 
particular important to understand the properties of this 
transition. 

In this paper we will calculate the gap around the 
limit of isolated rungs. We obtain reliable results in a 
wide range of parameters belonging to the rung-singlet 
phase. The transition curve to the spontaneously dimer- 
ized phase is computed. In addition, starting from the 
strong coupling limit of isolated rungs, the limit of iso- 
lated spin chains is discussed. 

We consider the S = i antiferromagnetic two-leg spin 
ladder plus additional four-spin exchange terms H cyc 
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The exchange couplings along the rungs and along the 
legs are denoted by J± and by Jy , respectively. J cyc de- 
notes the strength of the four-spin magnetic exchange 
terms. There is also another way based on cyclic permu- 
tations Pijki to include the leading four-spin exchange 
term. It differs in certain two-spin terms from Eq. (1) 

H P = JP ]T Si.lSi.2 + Jf\ S -- S < '•• + H cyc (2a) 
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Both Hamiltonians are identical except for couplings 
along the diagonals 3 if J± and J|| are suitably redefined. 
First, we use Hamiltonian H (l) since it is established 
that the four-spin terms are the most significant ones if 
the magnetic Hamiltonian is seen as effective model for 
the low-lying modes of a realistic insulating three-band 
Hubbard model 8 . But results for Hamiltonian H p (2) 
will also be presented. 

We use a continuous unitary transformation (CUT) 
to map the Hamiltonian H to an effective Hamiltonian 
H c g which conserves the number of rung-triplets, i.e. 
[H ,H eS ] = where H := i% ||=CJ , Jcyc=0 ] 18 - The ground 
state of iJ e ff is the rung-triplet vacuum. The effective 
Hamiltonian H c g is calculated in order 11 in x := J\\/ J± 



and x cyc := J C y C /Jx- Thereby, we obtained the ground- 
state energy Eq = (0\H c s\0) and the one-triplet disper- 
sion uj(k) — {k\H c fi\k) — Eq. The one-triplet dispersion 
ui(k) has a minimum for k = 7T, the one-triplet gap 
A (x, x cyc ) := u>(tt). By such perturbative approaches 
working on the operator level the spin ladder without 
cyclic exchange has been investigated previously with 
great success 19 ' 20 . 

The standard approach to calculate a phase transition 
line with series expansions is to use dlogPade approxi- 
mants on A(x, x cyc ). This yields reliable results only in 
a very small region about the exactly known phase tran- 
sition point [x — 1/5, x cyc = 1/5] (See grey square in 



Fig. 3 or similarly in Fig. 4). Generally, for x 
dispersion and the gap are known exactly 
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The results extrapolated in x are reliable for x € [0.1, 0.3] 
where the gap closes linearly in x and x cyc . For Hamil- 
tonian H p (2) the analogous situation is found at and 
about the exact point x p = x p yc = 1/4 as shown in Ref. 
21. Note that we use the parameters with superscript p 
to distinguish results for the Hamiltonian H p (2) clearly 
from those for the Hamiltonian H (1). 

In the following, we advance a recently introduced ex- 
trapolation technique 22 in order to investigate the rung- 
singlet phase for larger/lower values of x and x cyc . The 
main idea is to express the series expansion not in exter- 
nal parameters of the system like x and x cyc , but in an 
internal energy. Thereby, we combine high series expan- 
sion and renormalization group ideas. The natural inter- 
nal energy scale of the two-leg ladder is the one-triplet 
gap. In practice, we define the function 



G{x) = 1 - A(x) = 1 - 



A(x, rx) 
(1 + x) J± 
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where r = x cyc /x = J cyc /Ju will be kept constant for 
the extrapolation in x. The function G{x) behaves like 
G oc x for x — > so that any expansion in x can be 
converted in an expansion in G. Using the expansion 
for A(x) we calculated the inverse function x — x(G) as 
a series in G up to order 11 from Eq. (4). The quantity 
A = A/[(l+x)J±] measures the gap in units of J± + J\\ to 
ensure empirically a monotonic behavior of A as function 
of x. Then the existence of the inverse x{G) is assured. 
Next we consider the derivative of A(x) 
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dx 
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Substituting x = x(G) in Eq. (5) we obtain 



dG 
dx 



P(G) . 
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where P(G) is the truncated series of order 10 in G. Note 
that even the convergence of the truncated series P(G) is 
significantly better than the convergence of the truncated 
series A'(a;) in x as discussed in Ref. 22. Because the gap 
is a monotonic decreasing function for r =const we can 
use dlogPade approximants for P(G) since —dG/dx is 
non-negative. Integrating Eq. (6) yields 



dG 
PjGj 



dx 



x 



Therefore, integrating the left hand side to G = 1, i.e. 
A — 0, provides the phase transition point [xq, rxo] for a 
given r. For any Go S [0, 1[ the gap is A(xq, rx )/J± = 
(1 + xq)(1 — Go). In this way, A(x, x cyc ) is obtained. 

First, we examine the behavior of the gap in the limit 
of small r and G = 1. This corresponds to the situation 
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FIG. 1. For Hamiltonian H (1); black lines show the dlog- 
Pade approximants for dA/dx at G = 1 as a function of 
r = Xcyc/x. The grey line is a fitted spline which follows the 
asymptotic behavior Eq. (10) with A = 0.41 and A' = 0.85 at 
small values of r and approximates the available dlogPade re- 
sults. The points marked by stars are set by hand to guide the 
spline smoothly in the intermediate region. The extrapola- 
tions in Figs. 3 and 4 require actually only the values r < 0.5. 



of two spin chains which are weakly coupled by the four- 
spin interactions. Bosonization results show that the only 
relevant operator is the four-spin leg-leg interaction 21 . 
The triplet gap scales as 



XJ± — X Jq- 



(8) 



in leading order in J± and J cyc - Here A and A' are 

(5) 

non-universal constants 15 . In our case we have a criti- 
cal theory with central charge c = | and SU(2) symme- 
try which is described as the k — 2 Wess-Zumino-Wittcn 
model 15,23 . Rearranging Eq. (8) we obtain 



A 



X x c — x 
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where x c = A/(A'r) is the value of x where the gap van- 
ishes for given r. Therefore, the derivative of A for small 
r at G = 1, i.e. x = x c , is given by 



A(G = 1, 



(A'r) 2 
X + X'r 



(10) 



-A/x 2 



(7) In the case of r — > we expect A = and A 

— A /A from Eq. (9). Exploiting A = — A /A in a biased 
dlogPade approximant we find A = 0.4±0.03 in very good 
agreement with Quantum Monte Carlo results Aqmc — 
0.41 ± 0.01 from Ref. 24. 

In Fig. 1 the solid line corresponds to the dlogPade 
[7,2] for A (G = l,r). For r < 0.3 the asymptotic for- 
mula (10) is well reproduced by the approximant. A 
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FIG. 2. Same as in Fig. 1 for Hamiltonian H p (2) 

minute (not discernible) offset at r = occurs as a nat- 
ural consequence of the dlogPade approximation which 
describes a quantity of a given sign only. Using the value 
A = 0.41 ± 0.01 we deduce for the second non-universal 
constant A' the value 



A' = 0.85 ±0.2 



(11) 



If we perform the same analysis for Hamiltonian (2) we 
obtain Fig. 2 leading to the same result for A' given in 
Eq. (11). This fact corroborates the validity of the analy- 
sis and agrees perfectly with the finding by Miiller et al. 21 
stating that the relevant term in the cyclic exchange is 
the leg-leg coupling so that both Hamiltonians (1,2) lead 
to the same result for large leg couplings and small cyclic 
exchange couplings. 

For larger values of r or r p we interpolate between vari- 
ous approximants. This works better for Hamiltonian (2) 
(see Fig. 2) than for Hamiltonian (1) (see Fig. 1). But 
the interpolating functions are in any case quite similar. 
The uncertainty in the interpolation leads to the error 
bars in the subsequent extrapolations shown in Figs. 3 
and 4. These extrapolations are done for values r ;$ 0.5 
by subtracting the interpolated values depicted in Figs. 3 
and 4 from the truncated series for A (G) so that we ob- 
tain the series of a function that vanishes at G = 1. We 
find that many in this way biased dlogPade approximants 
yield reliable results. This supports our approach to in- 
clude the properties of the weakly coupled chains in the 
extrapolations. Finally the subtracted bias is re-added 
to arrive at the proper result. 

In the limit x — > oo, we conclude from Eq. (9) that the 
transition line converges against the asymptotic line 
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= A/A' re 0.52 ±0.14 
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FIG. 3. Extrapolated spin gaps for Hamiltonian -ff(l) in 
the [x, a; C y C ]-plane (see main text). The grey line is the ob- 
tained phase transition line A = and the grey square is the 
exactly known transition point [x — l/5,x cyc = 1/5]. The 
points marked by grey circles and error bars indicate the esti- 
mated accuracy of the extrapolations. On the left side of the 
transition line the system is in the rung-singlet phase, on the 
right side in the spontaneously dimerized phase. 



In Fig. 3 the extrapolated values of the spin gap of 
the Hamiltonian (1) in the [x, x cyc ]-plane are presented. 
The black solid lines denote A(xo, x cyc ) for a fixed x as 
a function of x cyc . These lines are shifted by xq in re- 
direction producing a quasi three-dimensional plot. The 
end-point of a black line corresponds to A(x,a; cyc ) = 0. 
These points yield the grey solid line which is the tran- 
sition line between the rung-singlet phase and the spon- 
taneously dimerized phase. As discussed above, we use 
biased approximants in the range x € [0.3, oo[ for the 
transition line. In the range x G [0.1,0.3] the unbiased 
extrapolations are safe due to the good convergence of the 
series obtained near the exactly known transition point 
(grey square). In the limit x — > even the truncated 
series gives quantitative results. Using Eq. (7) one finds 
in addition strong evidence for 
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using the values for A and A' obtained above. This re- 
sult holds again for both Hamiltonians (1,2). We cannot 

confirm the value of x^. mpt = 0.22 advocated in Ref. 21. ment with the results in Ref. 21. This is the reason why 



at x = where r) = 0.3 ± 0.02. The transition point, i.e. 
A = 0, for x = is found to be [0, 0.3 ± 0.002]. 

The smooth connection between the different extrapo- 
lations corroborates the reliability of our results in a wide 
region in the [x, x cyc ] -plane. 

In Fig. 4 the corresponding results for spin gap of the 
Hamiltonian 2 in the [x p , a;P yc ]-plane are depicted. The 
biased extrapolation is used for x p gl 0.4. Besides quan- 
titative differences there occurs one qualitative difference 
at low values of x p . For Hamiltonian (2) no closing of the 
gap on increasing x P yc for x p iS 0.1 was found in agree- 
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FIG. 4. Same as in Fig. 3 in the [x p , x%y C ] -plane for Hamil- 
tonian H p (2). The grey square is the exactly known transi- 
tion point [x p — 1/4, x P yc = 1/4]. The black circles are points 
taken from the curves in Ref. 21. The triangles are DMRG 
results (downward from Ref. 25; upward from Ref. 16). The 
diamond is determined from the maximization of the central 
charge by exact diagonalization 23 . 



the grey line is not prolonged below x p = 0.1. Apart 
from this point, the shape of the transition line is similar 
for both Hamiltonians. 

Quantitatively, it is interesting to compare to results 
obtained by other approaches, see Fig. 4. Refs. 25, 16 use 
density matrix renormalization. Another work analyzes 
the finite size scaling to determine the maximum central 
charge c. The cyclic exchange at which it appears deter- 
mines the critical x P yc , see Ref. 23. The spread and the 
error bars of the various results allow to assess the accu- 
racy of the data. We conclude that it is at present not 
yet settled where precisely the transition between gapped 
rung-singlet and dimerized phase occurs for the isotropic 
ladder. We propose to carry out careful finite-size scaling 
on data for periodic systems to clarify this issue. Note 
that the symmetry change between the rung-singlet and 
the dimerized phase cannot be represented properly in a 
single open system. 

In summary, we have investigated the rung-singlet 
phase of the S = ^ two-leg spin ladder with additional 
four-spin interactions. We used a continuous unitary 
transformation to calculate the one-triplet gap in a high 
order series expansion about the limit of isolated rungs. 
The use of an internal energy scale as the new expan- 
sion variable enabled us to calculate the transition line 
between the rung-singlet phase and the spontaneously 
dimerized phase reliably in a wide region of parameter 
space. Our results are consistent with the field theoretic 
results in the limit of weakly coupled chains. We repro- 
duce properties of the bosonization results in that limit. 
In addition, we give an estimate of the non-universal con- 



stants A and A' which appear in these bosonization treat- 
ments. The value for A is in very good agreement with 
quantum Monte Carlo results 24 . We have given an exam- 
ple that the combination of high order series expansion 
and renormalization group ideas can be a powerful tool. 

We thank A. Biihlcr, C. Knetter, U. Low and 
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